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ABSTRACT 

At the earliest epochs of structure formation in cold dark matter (CDM) cosmologies, the smallest 
nonlinear objects are the numerous small halos that condense with virial temperatures below ~ 10 4 K. 
Such "minihalos" are not yet resolved in large-scale three-dimensional cosmological simulations. Here we 
employ a semi-analytic method, combined with three-dimensional simulations of individual minihalos, 
to examine their importance during cosmological rcionization. We show that, depending on when reion- 
ization takes place, they potentially play an important role as sinks of ionizing radiation. If reionization 
occurs at sufficiently high redshifts (z r ^ 20), the intergalactic medium is heated to ~ 10 4 K and most 
minihalos never form. On the other hand, if z r <, 20, a significant fraction ( <; 10%) of all baryons have 
already collapsed into minihalos, and are subsequently removed from the halos by photoevaporation as 
the ionizing background flux builds up. We show that this process can require a significant budget of 
ionizing photons; exceeding the production by a straightforward extrapolation back in time of known 
quasar and galaxy populations by a factor of up to ~ 10 and ~ 3, respectively. 

subject headings: cosmology: theory - early universe - galaxies: evolution 



1. INTRODUCTION 

The lack of any Gunn-Peterson troughs in the spectra 
of distant quasars (Stern et al. 2000; Fan et al. 2000), 
as well as the presence of Lyman a emission lines in the 
spectra of high-redshift galaxies (Weymann et al. 1998; 
Hu et al. 1999), imply that the hydrogen in the intergalac- 
tic medium (IGM) is highly ionized by redshift z m 6. It is 
widely believed that reionization was caused by an early 
population of either galaxies or quasars. The process of 
reionization and its impact on several key cosmological is- 
sues have recently received much theoretical attention, and 
have been studied by several authors using semi-analytic 
models (Meiksin & Madau 1993; Shapiro et al. 1994; 
Haiman & Loeb 1997, 1998) and three-dimensional nu- 
merical simulations (Gnedin & Ostriker 1997). In general, 
these works aim to follow the time evolution of the fill- 
ing factor of ionized (HII) regions, based on some input 
prescriptions for the emissivity and spectra of the ionizing 
sources. 

More recent works have focused on the increased rate of 
recombinations in a clumpy medium relative to a homo- 
geneous one, i.e. when (p 2 ) > (p) 2 (Ciardi et al. 2000; 
Benson et al. 2000; Chiu & Ostriker 2000; Gnedin 2000a; 
Madau, Haardt, & Rees 1999 [hereafter MHR]). These 
studies have left significant uncertainties on the details of 
how reionization proceeds in an inhomogeneous medium. 
Since the ionizing sources are likely embedded in dense re- 
gions, one might expect that these dense regions are ion- 
ized first, before the radiation escapes to ionize the low- 
density IGM. Alternatively, most of the radiation might 
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escape from the local, dense regions along low column den- 
sity lines of sight. In this case, the underdense 'voids' are 
ionized first, with the ionization of the denser filaments 
and halos lagging behind (Miralda-Escude et al. 2000). 

In this paper, we point out that that the density inho- 
mogeneities at the earliest redshifts are dominated by the 
smallest nonlinear structures, i.e. halos near the cosmolog- 
ical Jeans mass, Mj eans w 10 4 M©. Such small scales have 
not yet been resolved in numerical simulations of reioniza- 
tion (e.g. Gnedin 2000a), and have also not yet been fully 
quantified in the semi-analytic works (although see Chiu 
& Ostriker 2000 and Benson et al. 2000 for partial treat- 
ments). Our main goals in this work are (1) to quantify 
the importance of the high-redshift minihalos as sinks of 
ionizing photons; and (2) to assess whether by extrapolat- 
ing to early times the known population of galaxies and 
quasars a sufficient number of UV photons are produced 
for hydrogen reionization. 

This paper is organized as follows. In § 2, we summarize 
the ionizing photon budget from the known populations 
of galaxies and quasars. In § 3, we describe our model of 
an individual photoevaporating minihalo, based on three- 
dimensional numerical simulations. Combining this model 
with a hierarchical structure formation scenario, in § 4 we 
explicitly show that the ensemble of minihalos dominate 
the clumping of the high-redshift IGM. In § 5, we then 
compute the total number of ionizing photons consumed 
by the photoevaporating minihalos in the same hierarchi- 
cal cosmology. In § 6, we argue that the minihalos must 
indeed have been photoevaporated. In § 7, we discuss the 



Minihalos and Reionization 



relevance of the covering factor of minihalos around ion- 
izing sources, and the uncertainties this implies for our 
results. Finally, in § 8, we offer our conclusions and sum- 
marize the implications of this work. Unless mentioned 
otherwise, throughout this paper we adopt a flat ACDM 
cosmology with (fi A , fi b , h, a 8 ,n) = (0.7, 0.04, 0.7, 0.9, 1). 

2. IONIZING PHOTONS FROM KNOWN SOURCES 

The two most natural candidates for reionization are ex- 
tensions of the known populations of quasars, or of galax- 
ies with ongoing star-formation, to rcdshifts beyond z > 6. 
Here we briefly review the number of ionizing photons ex- 
pected from these sources. This will serve as a reference 
point for comparison to the ionizing photon budget we 
obtain in our calculations below. 

2.1. Quasars 

The luminosity function (LF) of quasars has been mea- 
sured in optical surveys, and convenient parametric fit- 
ting formulae have been published by, e.g. Pei (1995). A 
more recent determination by the Sloan Digital Sky Sur- 
vey extends the measurement of the bright end of the LF 
to higher redshifts, and is consistent with earlier optical 
results (Fan et al. 2000). Here we adopt the fitting for- 
mulae of MHR, which has a somewhat shallower slope to- 
wards higher redshifts than, e.g., Pei (1995), and therefore 
predicts a larger number of quasars when extrapolated to 
high z. We have used this LF, together with the aver- 
age intrinsic spectrum of Elvis et al. (1995), to compute 
the production rate of ionizing (E > 13.6 eV) photons 
from quasars. The faint-end slope of the empirical LF is 
dlog<I>/<ilogZ/ = —1.64, so that most ionizing photons are 
produced by relatively bright quasars near the "knee" of 
the LF. A recent determination of the high-redshift LF 
in the soft X-ray band (Miyaji et al. 2000) has yielded 
significantly weaker evolution out to z « 4.5 than seen in 
the optical. Although the X-ray LF still has large uncer- 
tainties at z <; 2, taking the X-ray results at face value 
may indicate the existence of a larger abundance of high- 
redshift quasars than seen in the optical. If this inference 
is confirmed in future studies, it could increase the number 
of ionizing photons produced by bright quasars relative to 
our estimates below, unless this high-,? X-ray population 
is intrinsically obscured at optical/UV wavelengths. 

2.2. Galaxies 

The total star formation rate in the universe can be es- 
timated using the sample of high-redshift galaxies found 
by the Lyman-break technique (Steidel et al. 1999; Madau 
& Pozzetti 2000). Since the rest-frame UV continuum at 
1500-2800 A (redshifted into the visible band for distant 
sources) is dominated by the same short-lived, massive 
stars which are responsible for the emission of photons 
shortward of the Lyman edge, the needed conversion fac- 
tor, about one ionizing photon for every 5 photons at 1500 
A, is fairly insensitive to the assumed IMF and is indepen- 
dent of the galaxy history for t > 10 73 yr (MHR). We 
normalize the number of ionizing photons to the observed 
1500 A flux (rest- frame) , i.e. bypass the need for any cor- 
rection due to dust extinction (this is a good approxima- 
tion provided the mean color excess -E912-1500 = 1-6-Eb-v 
is actually small) . We have further assumed an average es- 
cape fraction of / csc = 50% for the ionizing radiation from 



the galaxy HI layers into the IGM, relative to the escape 
fraction at 1500 A. This value is a factor of ~ 5 higher 
than that inferred in nearby starbursts (Leitherer et al. 
1995). Theoretically, the escape fraction is expected to be 
lower at high redshifts, because of the higher gas densities 
then (Wood & Loeb 1999; Dove et al. 2000). However, 
a recent measurement of the mean escape fraction from a 
sample of 29 Lyman break galaxies at (z) = 3.4 has yielded 
a value > 50% (Steidel et al. 2000). Although the physi- 
cal implications of this result are still unclear (Haehnelt et 
al. 2000), to be conservative we adopted this value here. 
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Fig. 1. — Production rate (lower curves) and total number (inte- 
grated over cosmic time, starting from high redshift, upper curves) 
of H-ionizing photons from known quasars and galaxies, extrap- 
olated to z > 4.5. The calculations are based on the counts of 
high-redshift quasars and Lyman break galaxies. We assumed the 
escape fraction of ionizing photons to be / esc = 1 for quasars and 
/esc = 0.5 X / eS c, 1500 f° r galaxies. See text (§ 2) for discussion. 

The ionizing photon production rates for both quasars 
and galaxies are shown in Figure nl in units of photons 
above 13.6 eV emitted into the IGM per 10 8 years per in- 
tergalactic hydrogen atom. Also shown is the total num- 
ber of photons radiated per H-atom prior to redshift z. 
To obtain these integrals, we have extrapolated the photo 
emission rates dn^/dz beyond z = 6, using the observed 
slope dN/dz at z = 4. These extrapolations are clearly 
very uncertain. For reference, we note that if the photon 
emission rates had stayed flat at their z — 6 values be- 
yond z > 6, than in the stellar case, the total number of 
photons would increase very little (~ 20%), since the ob- 
served rate is already nearly constant in redshift. In the 
quasar case, this would result in an increase in the total 
number of photons generated prior to z — 6 by a factor 
of w 4. Under these extrapolations, as this figure reveals, 
luminous quasars produce little ionizing radiation prior to 
redshift z = 5, while galaxies provide « 4 x (/ csc /0.5) 
ionizing photons per hydrogen atom by this epoch. 

3. PHOTOEVAPORATION OF MINIHALOS 

In cold dark matter (CDM) cosmologies, the mean den- 
sity and temperature of the IGM imply the cosmological 
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where Mj is the total (gas + DM) mass of perturbations 
allowed to grow in linear theory. This mass corresponds 
to a virial temperature of the dark matter halo 
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K. 

It is natural to expect that the earliest nonlinear ob- 
jects are hosted in halos of the size described by equa- 
tions ([j]) and (H). Numerical simulations find that the 
mass scale can be up to a factor of ~ 10 still lower than 
this (Gnedin 2000b); however, as we shall find below our 
results are insensitive to the exact value of the smallest 
allowed halo size, as long as it is below T « 1000K. We 
shall hereafter refer to halos with virial temperatures be- 
tween Tj < T v i r < 10 4 K as "minihalos". A sufficiently 
strong early X-ray background (XRB) could catalyze the 
formation of H2 molecules in minihalos, which could then 
cool and fragment into stars, or form central black holes 
(Haiman, Rees & Loeb 1996; Haiman, Abel & Rees 2000). 
In this case, the gas from the minihalos would likely be 
ionized and expelled by the internal UV source. Alterna- 
tively, in the absence of any early XRB, most minihalos 
would likely not have sufficient molecular abundance to 
cool and dissipate (Haiman, Rees & Loeb 1997; Haiman, 
Abel & Rees 2000). In the further absence of any exter- 
nal UV radiation, they would then remain in approximate 
hydrostatic equilibrium until they eventually merge into 
more massive objects. 

However, photoionization by either a nearby external 
UV source, or by the smooth cosmic UV background ra- 
diation after the reionization epoch, heats the gas inside 
a minihalo to a temperature rs 10 4 K. By definition, the 
gas inside a minihalo is then no longer bound, leading to 
the photoevaporation of the baryons out of their host ha- 
los (Rees 1986; Shapiro, Raga & Mellema 1998; Barkana 
& Loeb 1999). In general, the resulting photoevapora- 
tive outflow has a complex velocity, density and ionization 
structure, whose time-evolution has been studied in a two- 
dimensional simulation by Shapiro & Raga (2000). In the 
present paper, we shall not address this evaporation pro- 
cess itself in detail. Our main goal here is to obtain a rough 
estimate of how many recombinations each hydrogen atom 
experiences during the evaporation process. We will there- 
fore use a simplified model of the evaporative flow, which 
we calibrate using three-dimensional simulations. 

3.1. Analytic Model of Photoevaporation 

We begin by adopting a radial density profile njj(r) for 
the undisturbed gas, based on the spherically symmet- 
ric truncated isothermal sphere (TIS, Shapiro, Ilyev & 
Raga 1999), with the scalings converted to our adopted 
ACDM cosmology as given in Navarro, Frenk & White 
(1997, hereafter NFW). This profile has the advantage of 
characterizing a spherical object with a well-defined core, 



calibrated specifically for use within the Press-Schechter 
formalism (which we shall use below). For this density 
profile, the average internal overdensity of the gas within 
the halo, relative to the background hydrogen density, is 
<5i nt = (nn)/n-R = 130.5. Similarly, the mean "internal 
clumping" is C[ nt = (i!-h)/"h = 444 2 . As an alterna- 
tive approach, we have assumed the gas to be in hydro- 
static equilibrium within a halo that has a density profile 
as parameterized by NFW, with a concentration parame- 
ter c = 5. This value was obtained by following the algo- 
rithm described in the appendix of NFW, which assigns 
to each halo of a given mass identified at redshift z a col- 
lapse redshift z co \\ (defined as the time at which half of 
the mass of the halo was first contained in progenitors 
more massive that some fixed fraction of the final mass), 
and assuming that the characteristic halo density is pro- 
portional to the critical density at z co n (Navarro, Frenk & 
White 1997). With this algorithm a minihalo of total mass 
M = 10 7 fc _1 M Q at z = 10 has (z co ii,c) = (13.5,4.8) (see 
also Madau, Ferrara, & Rees 2000). The NFW profile, by 
definition, has a mean internal overdensity of 200, and the 
embedded gas is somewhat more centrally condensed than 
the TIS we adopted, C; n t = 532 2 . We conclude that an 
NFW profile would give similar results to those obtained 
below, except with the total number of recombinations in- 
creased by about (532/444) 2 « 50%. 

Within the core of the TIS profile, the central density 
contrast reaches a value as high as ~ 2 x 10 4 relative to the 
cosmic background. In reality, for gas at temperature T V1I 
that has collapsed adiabatically, the density contrast can 
not exceed (T vir /Ti G M) 3 / 2 , where Tiqm is the IGM tem- 
perature (given by T lGM (z) w 2.73(l + z c )[(l + z)/(l + z c )] 2 
with z c w 150). We therefore modify the inner density pro- 
file, and adopt p(r) = min[p T i S (r), pi G M(T vir /Ti GM ) 3/2 ], 
where ptis(^) is the original TIS density run, and pigm 
is the density of the uniform IGM. This modification re- 
duces the average density contrast and the internal clump- 
ing below their fiducial values of 130.5 and 444 within ha- 
los whose virial temperature is lower than T v ; r <, 1500K. 
Below we find that w 68% of the recombinations in the 
ensemble of minihalos arise in those with T v ; r < 1500K, 
so that uncertainties in the density profiles of the smallest 
halos do not significantly effect our results. 

We next assume that the spherically symmetric halo is 
suddenly photoionized by an external UV flux. Simulta- 
neously, the gas temperature jumps to 10 4 K, resulting in 
an outflow. We assume that photoevaporation requires a 
sound crossing time, t pc ~ i? v ir/10 km s _1 , so that the 
photoionized gas retains its original shape for a duration 



tpe / 



Rv 



10 km s 



(3) 



after which it is fully dispersed into the IGM (where i? v i r 

is the radius of the TIS). Note that i pe ex T vi ( , so that 
smaller halos are photo-dissociated more rapidly. This 
prescription is admittedly oversimplified, since the gas 
should expand gradually, and non-uniformly, at varying 
speeds, rather than as a step-function. In addition, for 
low UV fluxes, the expansion might start before the gas is 
ionized (see discussion below). Nevertheless, as we will see 
below, the numerical simulations show a fairly rapid evap- 
oration, and can be used to calibrate the normalization 
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constant "/" that we have introduced. 

3.2. Numerical Simulations 

Our main goal is to estimate the number N pc of re- 
combinations each hydrogen atom experiences during the 
photoevaporation process. This number is given by 



iVpo — dt dVx 2 n^aB X / dVnn 

= f -,m 7Cint"-HaB/^int- 

10 km s i 



(4) 



Here the volume integrals are taken over the region that 
contains the virial gas mass, x is the ionized fraction 
(which is a function of radius and time), and ub = 
2.6 x 10~ 13 cm 3 s _1 is the case B hydrogen recombination 
coefficient evaluated at T = 10 4 K. Given an initial density 
profile, the second line can be readily computed, except 
for the normalization constant /. For reference, we note 
that the photoevaporation time (eq. ||) for a T v „ — 10 4 K 
minihalo at redshift z — 10 is approximately 25/% of the 
Hubble time at that redshift; we find that during this in- 
terval each hydrogen atom within this halo recombines 
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340/ times. 



The first line in equation (3) is computed explicitly in 
the simulation, which gives the total number of recombi- 
nations each hydrogen atom experienced during the whole 
photoevaporation process. By comparing this number to 
the second line, we determine the normalization constant 
/. The simulations are performed by setting up stable 
truncated isothermal spheres in both dark matter (DM) 
and gas as the initial condition for the three-dimensional 
structured adaptive mesh refinement cosmological hydro- 
dynamics code enzo of Bryan & Norman (1997; 1999). For 
simplicity, we use a standard CDM cosmology with h = 0.5 
and Ob = 0.1. We solve the chemical reactions for the 
ions of hydrogen and helium, including collisional ioniza- 
tion, radiative recombinations, and photo-ionization. The 
reaction rates are taken from Abel et al. (1997). The hy- 
drogen photo-ionization rate is assumed to be uniform on 
the grid(s); i.e. we do not solve radiative transfer. The 
photo-ionization heating rate is simply given by fc p h(e p h), 
where (e p h) denotes the average energy of photo-electrons 
(which depends on the ionizing spectrum; here we adopt 
<e ph ) = 2eV). 

The typical photoionization rates experienced by mini- 
halos depend strongly on the luminosity and clustering 
properties of the radiation sources. In addition, the pho- 
toionization rates can vary substantially for individual 
minihalos, because of the 1/r 2 dependence of the flux from 
nearby sources. In our simulations, we have experimented 
with photoionization rates of fc p h = 10~ 13 — 10~ 10 s _1 , cor- 
responding roughly to a background flux of J = 3 x 10~ 23 — 
10~ 20 erg s _1 cm~ 2 Hz -1 sr _1 at 13.6eV. In general, we 
find that for fast photoionization rates (A; p h Jj 10~ 12 s _1 ), 
the gas is fully ionized before any motion occurs. In this 
regime, we find a simple scaling of the numerical results 
(see below) with halo mass and redshift, independent of 
the actual value of fc p h- This is because the heating- 
induced motions are independent of the photo-heating 
rate, as long as this rate is comparable to or shorter than 
the sound crossing time. In the following, we will quote 



results from simulations with a high-enough flux that al- 
lows us to use these scalings, and thus eliminate the need 
for a large grid of simulations in mass, redshift, and flux. 

We emphasize that for smaller UV fluxes, the num- 
ber of recombinations can be reduced relative to our 
adopted estimate. The core of the halo has a lower ion- 
ization fraction, both because the equilibrium fraction is 
reduced, and because the photoionization time scale in the 
core becomes comparable to the sound crossing time (for 
fcph ^ 10~ 12 s _1 ), so that equilibrium is not established 
before bulk motion occurs. In addition, radiative trans- 
fer effects for a low flux might keep the central regions 
shielded until the density is significantly reduced, and fur- 
ther reduce the number of recombinations. In general, full 
radiative transfer would be needed to adequately quantify 
these effects, which we defer to a subsequent paper. Here 
we only estimate these effects in the optically thin case, 
and find that for the range of /c p h of interest, the reduction 
in the number of recombinations is small. As an example, 
for the low rate of /c p h = 10~ 13 s _1 , the core of our fiducial 
z = 10 halo would reach x — 0.4, resulting in an overall re- 
duction of the internal clumping of the ionized gas within 
the halo by a factor of two. Finally, we note that prior 
to reionization, it is conceivable that, due to the strong 
fluctuations of the UV flux as a function of position, the 
gas condenses and photo-evaporates from DM minihalos 
multiple times, and hence absorb many more photons. 

Our calibrations are based on four simulations, using 
halo masses of M = 10 6 M Q or M = 5 X 10 6 M Q ; and 
initial redshifts of z = 10 or z = 20. In all of these 
cases, by our choice of the flux, photoevaporation occurs 
very rapidly, with most recombinations occurring within 
a sound-crossing time; justifying the assumption of our 
simple model in equation (H). As an example, in Figure 
we show the evolution of the M = 10 6 M© halo starting 
at z = 10. The figure reveals that by z = 9.5 the density 
contrast is everywhere reduced below <5 ^ 10. In each of 
the cases we examined, we infer / = 1 from equation (H) 
to within ~ 15%. For the M = 10 6 M halo, we find 
/ = 1.03,0.86 for z = 20 and z = 10, respectively. Sim- 
ilarly, for M = 5 X 10 6 M , we find / = 0.85,0.89 for 
z — 20 and z — 10. Generally, considerable differences 
in the inferred values of / can be expected if one were to 
carry out radiative transfer calculations for more realistic 
density distributions and source evolutions. In the rest of 
this paper, we take our four simulations at face value, and 
set / = 1 in the remaining calculations. 
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Fig. 2. — The evolution of the density profile of a halo with mass 
10 6 Mq, illuminated by a background flux at z = 10. 



4. MINIHALOS AND GAS CLUMPING 

In this section, we illustrate the importance of minihalos 
by estimating their contribution to the overall dumpiness 
of gas in the universe. For simplicity, we assume sudden 
reionization at some redshift z r , corresponding to the over- 
lap of discrete HII regions, and a sudden increase in the 
background flux at photon energy E = 13.6eV from zero 



to a value J = J21 x 10 21 erg s 1 cm 2 Hz 



sr 



Prior 



to z r the background flux is assumed to be negligible, so 
that minihalos form down to the Jeans mass, or down to 
the virial temperature of Xj (eq. 0). 

At the overlap epoch, the minihalos with T v ; r < 10 4 K 
contain a mass fraction 



Im{zt) = ■F co u(z r ,Tj) - F co n(z r , 10 4 K) 



(5) 



of all baryons, where F co \\(z, T) is the Press-Schechter col- 
lapsed fraction at redshift z above the mass cutoff cor- 
responding to the virial temperature T. The minihalos 
therefore fill a fraction 



fv(z r ) 



Jm Or) 

130.5 



(6) 



of the volume, where (5; n t = 130.5 is the average over- 
density within a minihalo relative to the background. The 
fraction we use in our calculations is somewhat larger than 
this value because, as discussed in § 3.1. above, we reduce 
the central densities of the smallest minihalos relative to 
the TIS model; as a result, they occupy a larger volume. 

Once the ionizing flux turns on, the minihalos are pho- 
toionized, and each minihalo contributes to the gas clump- 
ing and increases the universal mean recombination rate, 
until it is photoevaporated. To illustrate the importance 
of minihalos, it is useful to define an average gas clumping 
factor, 

K) 



where the three terms represent the contributions from the 
smooth IGM, halos with T v ; r > 10 4 K, and minihalos, re- 
spectively. Previous works have computed gas clumping 
in a similar manner, but have not included the contribu- 
tion C< from minihalos (see Benson et al. 2000). Each 
term in equation (Q) can be written as C{ — /v,i(^f)/wjj, 
where fv.i is the volume filling factor and m is the num- 
ber density of ionized hydrogen in the i th component. The 
first term is obtained by assuming that the uncollapsed gas 
mass fraction « [1 — F co \i(z,Tj)] is fully ionized, and uni- 
formly fills the available volume « [1 — F co n(z, Tj)/130.5] 
(in our calculations, we have included the small correction 
to these expressions due to the reduction of the central 
densities of small minihalos). Note that the uncollapsed 
gas fraction is underdense relative to the overall mean of 
the universe, so that Ciqm < 1- 




C 



ClGM + C> + C< 



(7) 



Fig. 3. — Average clumping factor of ionized gas as a function 
of reionization redshift, when minihalos are included (solid curve), 
and when only halos with T V \ Y > 10 4 K are included (dotted curve). 
The dashed curve is from a hydrodynamical simulation by Gnedin 
& Ostriker (1997). 

The contribution from large halos, C> is computed by 
summing over all halos with T v - lr > 10 4 K, whose abun- 
dance is assumed to follow the Press-Schechter mass func- 
tion. These halos are collisionally ionized when they viri- 
alize, and subsequently cool via atomic line cooling. In 
the absence of an external UV field, most of the hydro- 
gen would recombine, and remain neutral. When exposed 
to the photoionizing background, the large halos are self- 
shielding, and are still only partially ionized from the sur- 
face inward. This renders estimates of the contribution 
from large halos to the overall clumping highly uncertain, 
since the gas in each halo can continue to cool, self-shield, 
and collapse to form a disk, while it is being illuminated 
by the external radiation field. Nevertheless, in order to 
derive an illustrative estimate, we have computed the ra- 
dius n to which the halos are ionized from the outside in, 
assuming that the gas density profiles follow that of the 
TIS solution. We have assumed further that the spherical 
halo is exposed to an external ionizing flux of J21 = 1, 
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and repeated the ID radiative transfer calculation for sev- 
eral halo masses, to obtain fj as a function of Mh a i - We 
then evaluated (nf ) for each halo mass, excluding the neu- 
tral core < r < r^ in each case, and obtained the final 
clumping factor by averaging (nf) over Mhaio- 

Finally, the contribution C< from minihalos is computed 
by summing the internal clumping C- m t < AAA 2 over all 
minihalos with 7j < T v ; r < 10 4 K, assuming that each 
minihalo is fully ionized, but only exists for a time t po 
given by equation (||). To illustrate the importance of 
minihalos, the resulting clumping factors Cigm + C> + C< 
and Cigm+C> are shown by the solid and dotted curves in 
Figure 0, respectively. Also shown for comparison is the 
clumping factor in a similar ACDM cosmology from the 
numerical simulation of Gnedin & Ostriker (1997, dashed 
curve). Our calculations show C as a function of the as- 
sumed reionization redshift, while the simulation result is 
adopted from a specific reionization history; this simple 
comparison is therefore somewhat inappropriate. Never- 
theless, as the figure shows, the 3D simulation results are 
well reproduced by the ensemble of partially ionized spher- 
ical halos with T vlr > 10 4 K, assuming a fixed background 
flux. However, when the minihalos are included, clumping 
sets in at a significantly higher redshift (z ~ 20 instead 
of z « 10). Although the relative importance of miniha- 
los decreases at lower redshifts, when larger scales start to 
collapse, the clumping factor is still enhanced by a factor 
of ~ 2 at z k, 6. 

5. MINIHALOS AS PHOTON SINKS 

We will now explicitly compute the total number of 
ionizing photons consumed inside minihalos following the 
overlap epoch z r . This equals the number of recombina- 
tions that occur while the gas is evaporated out of the 
minihalos. Although larger halos with T v ; r > 10 4 K even- 
tually dominate the gas clumping, as emphasized in the 
previous section, large halos can self-shield, cool, and col- 
lapse to form disks, making it difficult to estimate their 
overall contribution. In order to avoid these complications, 
below we focus on minihalos, and simply ignore halos with 
T vir > 10 4 K. 

Consider a post-overlap redshift z < z r , with At(z, z r ) 
denoting the time elapsed between z r and z. By this red- 
shift, the smallest minihalos have been evaporated, but 
those with virial temperatures above 



T pc (z) = 10 4 



At(z, Z r ) 



V(^,10 4 K) 



K 



(8) 



are still present (cf. eq. ||) and are contributing to the gas 
clumping. The contribution of minihalos to the clumping 
factor at redshift z is accordingly given by 



C<(z) 



Af(10 4 K) 



dM 



M(T min ) 



f dN 

\dM 



M 



<5intPbg 



C i r 



(9) 



where the first term in the integrand is the Press-Schechter 
mass function; the second term represents the volume- 
filling factor of halos of mass M; and the last term is the 
internal gas clumping of a halo of mass M, Ci n t (M) < 444. 
The lower limit of the integration corresponds to the small- 
est halo that has not yet been photoevaporated, whose 
virial temperature is 



T mi n(z) = max(Tj,T pe ). 



(10) 



Since Tj is always evaluated at redshift z r , the only z— 
dependence of equation (M) are through the minimum tem- 
perature of minihalos still being evaporated, and through 
the decrease of the background density pb g - Some of the 
minihalos are destroyed by mergers; however, we find from 
the extended Press-Schechter formalism (Lacey & Cole 
1993; eq. 28) that in the relevant ranges of redshifts (cor- 
responding to the photoevaporation time) and halo masses 
this fraction is always small. The total number of recom- 
binations per hydrogen atom in the universe that takes 
places in minihalos between z r and z is then given by 



t(z) 

N ICC = I dta B C < {z)n Hfi (l + z) 3 , (11) 

lt(z r ) 



where uh,o is the present-day average density of hydrogen 
atoms. Figure shows the evolution of A roc , assuming dif- 
ferent values for the reionization redshift, 6 < z r < 16. In 
general, the number of recombinations associated with the 
photoevaporation of minihalos is significant. For popular 
values for the reionization redshift, the minihalos consume 
over 10 ionizing photons per hydrogen atom. Most recom- 
binations take place within a short time after reionization 
- once the minihalos are evaporated, they no longer con- 
tribute to recombinations, as demonstrated by the rapid 
flattening of the curves to their asymptotic values in Fig- 
ure ||. We recall that smaller halos are photoevaporated 
more rapidly than larger ones, which makes their impor- 
tance diminish towards lower masses, despite their increas- 
ing abundance. We find that « 68% of the recombinations 
shown in Figure [| take place in those with T v ; r <S 1500K; 
and excluding halos with T v ; r < 1000K would make a neg- 
ligible difference. Our results imply that if the minihalos 
were indeed photoevaporated (see discussion below), the 
ionizing sources must have produced at least ~ 10 ion- 
izing photons per background hydrogen atom. This is a 
necessary, but not a sufficient condition for reionization 
(e.g. a sudden burst of 10 ionizing photons per H atom 
could evaporate the minihalos and reionize the universe, 
but would not ensure that most of the volume is subse- 
quently kept ionized) . 
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the existing minihalos disappear by merging, as they be- 
come parts of larger halos. Assuming the minihalos to be 
neutral and randomly distributed in space, the comoving 
photon mean free path d at redshift z is given by 



1 + z f M ° 

—i-= dMn(M,z)a(M)f smv (M,z r ,M ,z) 

a J Mj 

rM 

dM , 
im 3 \dM, 



dN 



(l + z) 3 (7Ti?2. r )/ surv , (12) 



where a(M) is the geometrical cross-section, i? v j r is the 
virial radius, n(M, z) is the physical abundance at redshift 
z, and (dN/dM) Zr is the comoving abundance at redshift 
z r of minihalos of mass M; and f SU rv(M, z r , Mq, z) is the 
probability that a halo of mass M at redshift z r has not 
become part of a halo of mass > Mq by redshift z. Here 
and in the upper limit of the integrals in equation ( |l2| ) , Mq 
is chosen to correspond to a virial temperature of 10 4 K. 
The survival fraction is computed as 



Fig. 4. — Total number of recombinations inside minihalos, follow- 
ing sudden reionization at redshift z r . Recombinations are shown 
per hydrogen atom, assuming the mean IGM density (see eq. |ll| ). 

An interesting feature in Figure is the dependence of 
the number of recombinations on the assumed reioniza- 
tion redshift. If reionization occurs early (z r <; 20), i.e. 
by efficient ionizing photon production in rare, high-tr ob- 
jects, the abundance of minihalos is then still low. Since 
we do not allow minihalos to condense after reionization, 
most minihalos never form. As a result, minihalos arc less 
important sinks of ionizing radiation. Likewise, if reioniza- 
tion occurs late (closer to z — 6), then the minihalos have 
started merging into larger halos with T v ; r > 10 4 K. In this 
case, a fraction of minihalos reside in deep potential wells 
and can be self-shielded, reducing their role as photon 
sinks. The effect of minihalos are maximal if reionization 
occurs around the popular values of z rs 8 — 10, the red- 
shift when the minihalo abundance also peaks (note that 
this redshift depends sensitively on the assumed cosmol- 
ogy, and on the normalization of the power spectrum a&). 

In our treatment, we assumed sudden reionization. The 
later stages of reionization, corresponding to the slow 
outside-in ionization of dense regions on larger scales, can 
last a considerable fraction of the Hubble time. However, 
the initial overlap of HII regions likely proceeds rapidly 
(Haiman & Loeb 1998; Miralda-Escude et al. 1999; MHR; 
Gnedin 2000a). Nevertheless, one could follow the volume 
filling factor of ionized regions, and allow the formation 
of new minihalos only outside these regions (as is done 
in Haiman & Loeb 1998). Here we simply note that if 
the overlap epoch lasts between redshifts z\ and Z2, then 
reading N rcc from Fig [| at these two redshifts would ap- 
proximately bracket the number of photons consumed. 

6. THE MEAN FREE PATH OF IONIZING PHOTONS 

The minihalos that formed prior to the reionization red- 
shift z r contribute to the Lyman-continuum opacity of the 
universe and reduce the mean free path of ionizing pho- 
tons at later redshifts z < z r . No new minihalos form 
after reionization, and furthermore, a fraction 1 — / surv of 



fsurv(M,Z r ,M ,z) 



dMj_p(M, z r ,M b 



(13) 



M, 



where p(M, z r , Mi , z) is the probability that a halo of mass 
M at redshift z r is part of a halo of mass Mi ± dM\/2 at 
redshift z, given by Lacey & Cole (1993, eq. 28). 

We find that the mean free path is small compared to the 
Hubble distance ctHub(z)- For example, taking z r — 10 as 
the reionization redshift, we find d = 2.1 (comoving) Mpc 
at Z = 10. Equation (12) also shows that the contribution 
to the opacity from halos of mass M scales approximately 
as oc M -1 / 3 , implying that once the low-density IGM is 
ionized, but the halos are still neutral, the minihalos are 
the dominant source of opacity. As a result, the typical 
fate of an ionizing photon emitted around redshift z w 10 
is to be absorbed by a minihalo within a small fraction of 
the Hubble distance. 

Although most ionizing photons are absorbed by mini- 
halos, it is interesting to follow the evolution of the mean 
free path at z < z r under the assumption that the miniha- 
los stay neutral. This could be the case, e.g. if only ~ one 
ionizing photon is produced in the universe per H atom 
after z r . We find that approximately half of the minihalos 
present at z = 10 survive merging into large halos un- 
til z = 3, and the surviving minihalos lead to mean free 
paths of d — 10, 16, and 32 Mpc at z = 5, 4, and 3, re- 
spectively. This implies the existence of dN/dz » 54,42, 
and 31 minihalos per unit redshift at these redshifts. The 
minihalos considered here have hydrogen column densities 
of JVh *i 10 18 cm -2 , and their abundances should therefore 
be compared to that of high-redshift Lyman limit systems 
(LLS). Although the latter is not well measured at z « 5, 
it is inferred from high-resolution quasar spectra to be 
dN/dz » 2 ± 0.5 at z = 3 and dN/dz w 3.5 ± 1 at z = 4 
(e.g. Storrie-Lombardi et al. 1994), at least a factor of 
~ 10 lower than the number of surviving minihalos. We 
conclude that the observed relatively low number of LLS 
gives a strong constraint on the number of minihalos that 
could have survived photoevaporation. Unless future ob- 
servations reveal a 10-fold increase in the abundance of 
Lyman limit systems from z — 4 to z — 5, this provides 
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strong evidence that most minihalos were indeed photo- 
evaporated by z rs 5 (or at least their total geometrical 
cross-sections were reduced by an order of magnitude). 
Indeed, Abel & Mo (1998) have shown that halos with 
T v ; r <; 10 4 K already have a sufficient abundance to ac- 
count for all LLS at z <, 4; leaving no room for additional 
minihalos. 

7. THE COVERING FACTOR OF MINIHALOS 

So far we have assumed that reionization occurs sud- 
denly. Stated more precisely, in this picture reionization 
occurs in three stages. First, the HII regions around in- 
dividual ionizing sources rapidly overlap and fill most of 
the volume, occupied by the low-density IGM. Second, the 
minihalos are photoevaporated by the background radia- 
tion. Finally, the larger-scale dense regions are gradually 
ionized from the outside inwards. This scenario is justi- 
fied as long as the minihalos have a small covering factor 
F cov around an ionizing source. Here F cov is defined as 
the fraction of the 4n solid angle around a typical ionizing 
source occupied by minihalos out to a distance r s , where 
r s is the typical separation of the ionizing sources. On the 
other hand, if _F cov <; 1, then the minihalos must be photo- 
evaporated before the low-density regions can be ionized, 
since ionizing photons are necessarily absorbed in mini- 
halos before they are able to reach most of the volume. 
In this case, the first two stages of reionization would be 
reversed, and a typical minihalo is evaporated by a neigh- 
boring UV source, rather than by a uniform background. 
The resulting "lopsided" photoevaporation is considerably 
more complicated to model than the evaporation by a uni- 
form background, as we assumed in § 3. Although such 
modeling is beyond the scope of the present paper, we shall 
estimate the covering factor below. 

To estimate the covering factor, let us assume that the 
ionizing sources reside in halos with virial temperatures 
T vil rs 10 4 K, or mass Mq. Taking z r = 10, the average 
separation of these sources at the reionization epoch is 
r s = [MoidN/dMo)}- 1 / 3 = 0.5 (comoving) Mpc. If the 
minihalos were randomly distributed, F cov could be com- 
puted analogously to the mean free path (cf . eq. |lj) , and 
would be given simply by 



^cov,o - ^ ~ 23 ~ 25%. 



(14) 



However, equation (J14J) ignores the fact that miniha- 
los are strongly clustered around the ionizing sources, and 
therefore cover a larger fraction of the sky. Taking clus- 
tering into account, the abundance of minihalos of mass 
Mi at a distance r ±dr away from an ionizing source in a 
halo of mass Mq is given approximately by 



Namil + bobiUtr)], 



(15) 



where N% = Mi(dN/dMi) is the average minihalo abun- 
dance, £ m (r) is the mass correlation function at redshift 
z = 10, and b\ = b(z r ,Mi) and bo = b(z r ,M ) are the 
bias of halos of mass Mi and Mq relative to mass. The 
smallest spatial scale we consider is the virial radius of a 
halo at the Jeans mass, rs 1 (comoving) kpc, which is in 
the mildly non-linear regime. For simplicity, we use the 
linear correlation function and halo bias, although nonlin- 
ear effects could increase the clustering somewhat above 



the estimates below. The relevant halo masses are always 
above M* , i.e. even the smallest minihalo has a positive 
bias parameter b(z, M). 

For halos of a given mass Mi , clustering then increases 
the covering factor to 



^cov — -^ cov 



,o[l 



bpbi 



dr£ m (r)). 



(16) 



Rvi. 



Here i? v ; r refers to the virial radius of the minihalo of mass 
Mi . Averaging the term in square brackets over all masses 
between Mj and Mq, we find that the overall increase at 
redshift z = 10 is a factor of rs 3. The mass-dependence 
is primarily through the bias parameter b\, since the lin- 
ear mass correlation function £ m (r) is quite flat, and the 
radial integral is dominated by large radii. 

In summary, our estimates imply a typical covering fac- 
tor of F cov rs 75%. This relatively large value demon- 
strates that most ionizing photons emitted by a source 
are absorbed in minihalos before most of the volume of 
the IGM can be ionized. Our estimates are inherently 
uncertain, as they are based on spherical halos and the 
linear correlation functions. As higher resolution three- 
dimensional simulations become feasible in the future, it 
will possible to directly compute the covering factor of 
minihalos, and eventually to model their asymmetric pho- 
toevaporation (as done in two dimensions by Shapiro & 
Raga 2000). 

8. CONCLUSIONS 

We have shown that minihalos, i.e. halos with virial 
temperatures between the cosmological Jeans tempera- 
ture, Tj, and 10 4 K, play an important role during cosmo- 
logical reionization. In the earliest stages of reionization, 
the minihalos dominate the average gas clumping in the 
universe, and can therefore dominate the overall recom- 
bination rate. Depending on the reionization redshift, up 
to rs 20% of all baryons condense into minihalos before 
their formation is suppressed by the reheating feedback. 
The fate of the chemically pristine minihalos depends 
on the abundance of H2 molecules (see Haiman, Abel & 
Rees 2000). A sufficiently strong early X-ray background 
(XRB) could catalyze the formation of H2 molecules, and 
the minihalos could then cool and harbor internal ionizing 
sources, likely photo-evaporating their gas content from 
the inside (note that X-rays would have several additional 
important consequences for reionization, see Oh 2000). In 
absence of any early XRB, however, most minihalos would 
likely not have sufficient molecular abundance, and they 
would then remain in approximate hydrostatic equilibrium 
until photoevaporated by external ionizing sources. 

We find that the photoevaporation of existing minihalos 
in the range of virial temperatures 1000 K <, T V i r <, 10 4 K 
dominates the necessary ionizing photon budget, and re- 
quires 10 — 20 ionizing photons per background hydrogen 
atom. This value exceeds by about an order of magnitude 
the number of ionizing photons produced by an extrapo- 
lation of known populations of quasars to higher redshifts. 
An extrapolation of the Lyman break galaxy population to 
high redshift comes closer to meeting the required photon 
budget, provided that a large fraction of the ionizing pho- 
tons produced in these galaxies leak into the IGM. If the 
mean escape fraction is ~ 50% (cf. Steidel et al. 2000), 
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then the discrepancy is only a factor of ~ 3. Our con- 
clusions depend further on the type of sources responsible 
for reionization. Our results show that if the reionizing 
sources have typical separations of halos with virial tem- 
peratures of 10 4 K (or larger), then most minihalos are 
photoevaporated before most of the volume of the universe 
is reionized, because their covering factor around the ion- 
izing sources is of order unity. If the reionizing sources 
are more closely packed, in principle they could ionize the 
low-density regions, i.e. most of the volume, before the 
minihalos are photoevaporated. In either case, we find 
that the minihalos must have been photoevaporated by 
z w 5, in order not to overpredict the number of Lyman 
limit systems by a factor of ~ 10. 

In a previous study of inhomogeneous reionization, 
Miralda-Escude et al. (2000) have suggested that ~ one 
ionizing photon per hydrogen atom is sufficient to reionizc 
most of the volume of the IGM by z w 5. The gas clumping 
in that study was adjusted to match the results of numer- 
ical simulations, which did not resolve the smallest scales, 
i.e. scales on which minihalos dominate the gas clumping. 
Our results suggest that the necessary photoevaporation 
of dense gas out of minihalos can raise the photon require- 
ment at z Jj 5 by up to an order of magnitude. In arriv- 
ing at this conclusion, we have relied on a combination of 
a semi-analytic approach and three-dimensional simula- 
tions of individual halos, albeit without radiative transfer, 
and we have also made a number of simplifying assump- 



tions. More realistic cosmological hydrodynamical simula- 
tions that include realistic sources and three-dimensional 
radiative transfer will be available in the near future, and 
will be able to test our simplified model. 

The additional ionizing photons required for the photo- 
evaporation of minihalos could arise from an early popu- 
lation of low-luminosity quasars ("miniquasars"), whose 
abundance declines significantly less than that of optically 
bright QSOs. Alternatively, ionizing photons could have 
been produced in "minigalaxies" associated with low-mass 
halos collapsing at high redshift. Although both types of 
sources could have escaped detection with present instru- 
ments, both are well within the capabilities of direct imag- 
ing with future telescopes, such as the Next Generation 
Space Telescope. 

We thank Houjun Mo for useful discussions, and Jordi 
Miralda-Escude, Martin Haehnelt, Martin Rees, Rennan 
Barkana, and the referee, David Weinberg, for useful com- 
ments on the manuscript. We also thank Greg Bryan and 
Mike Norman for permission to use their adaptive mesh 
refinement code enzo. This research was supported by the 
NSF under Grant No. PHY94-07194 at the ITP, and by 
NASA through the Hubble Fellowship grant HF-01119.01- 
99A, awarded to ZH by the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities 
for Research in Astronomy, Inc., for NASA under con- 
tract NAS 5-26555. PM acknowledges support by NASA 
through ATP grant NAG5-4236. 



REFERENCES 



Abel, T., Anninos, P., Zhang, Y., & Norman, M. L. 1997, New Ast, 

2, 181 
Abel, T., & Mo, H. J. 1998, ApJ, 494, L151 
Barkana, R., & Loeb, A. 1999, ApJ, 523, 54 



Benson, A. J., Nusser, A , Sugiyam 



MNRAS, in press, astro-ph/0002457 
Bryan, G. L., & Norman, M. L. iy9V, in 



N., & Lacey, C. G. 2000, 



l Computational Astrophysics, 

ASP Conf. Ser. No. 123, eds. D.A. Clarke and M. Fall 
Bryan, G. L., & Norman, M. L. 1999, in Workshop on Structured 

Adaptive Mesh Refinement Techniques, ed. N. Chrisosoides, p. 167 
Chiu, W. A., & Ostriker, J. P. 2000, ApJ, 534, 507 
Ciardi, B., Ferrara, A., Governato, F., & Jenkins, A. 2000, MNRAS, 

314, 611 
Dove, J. B., Shull, J. M., & Ferrara, A. 2000, ApJ, 531, 846 
Elvis, M., Wilkes, B. J., McDowell, J. C, Green, R. F., Bechtold, J., 

Willner, S. P., Oey, M. S., Polomski, E., k. Cutri, R. 1994, ApJS, 

95, 1 . , 

Fan, X. etal. 2000, AJ, in press, astro-ph/000541^ 

Gnedin, N. Y. 2000a, ApJ, 535, 5311 

Gncdin, N. Y. 2000b, ApJ, 542, 535 

Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 

Haehnelt, M.,G,, Madau, P., Kudritzki, R., & Haardt, F. 2000,ApJ, 



submitted, 



astro-ph/0010631 
Haiman, Z., Abel, T., & Rces, M. J. 2000, ApJ, 534, 11 
Haiman, Z., & Loeb, A. 1997, ApJ, 483, 21 
Haiman, Z., & Loeb, A. 1998, ApJ, 503, 505 
Haiman, Z., Rees, M. J., & Loeb, A. 1996, ApJ, 467, 522 
Haiman, Z., Rees, M. J., & Loeb, A. 1997, ApJ, 476, 458 
Hu, E. M., McMahon, R. G., & Cowie, L. L. 1999, ApJ, 522, L9 
Lacey, C, & Cole, S. 1993, MNRAS, 262, 627 
Leitherer, C, Ferguson, H. C, Heckman, T. M., & Lowenthal, J. D. 

1995, ApJ, 454, L19 . 

MadarL EL, Ferrara, A., & Rees 2000, ApJ, submitted, astro 



ph/C 010158 



Madau, P., Haardt, F., & Rees, M. J. 1999, ApJ, 514, 648 [MHR] 

Madau, P., & Pozzetti, L. 2000, MNRAS, 312, 9 

Meiksin, A., & Madau, P. 1993, ApJ, 412, 34 

Miralda-Escude, J., Haehnelt, M., & Rees, M. J. 2000, ApJ, 530, 1 

Miyaji, T., Hasinger, C, & Schmidt, M. 2000, A&A, 353, 25 

Mo, H.J., Mao, S. & White, S.D.M. 1998, MNRAS, 295, 319 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 490, 493 

(NFW) . . 

Oh, S. P. 2000, ApJ, submitted, astro-ph/0005262 

Pei, Y. C. 1995, ApJ, 438, 623 ' 

Press, W. H., & Schechter, P. L. 1974, ApJ, 181, 425 
Rees, M. J. 1986, MNRAS, 218, 25p 

Shapiro, P. R., Giroux, M. L., & Babul, A. 1994, ApJ, 427, 25 
Shapiro, P. R., Ilyev, I., & Raga, A. C. 1999, MNRAS, 307, 203 
Shapiro, P. R,., & Raga, A. C. 2000, in Cosmic Evolution and 
Galaxy Formation: Structure, Interactions, and Feedback, ASP 
Conference Series, ^ds , J , Franco, E , .Terlevich, O. Lopez-Cruz, I. 



Aretxaga, in press, astro-ph/0004413 
Shapiro, P. R., Raga, A. C, & Mellema, G. 1998, in Molecular 

Hydrogen in the Early Universe, Memorie Delia Societa 

Astronomica Italiana, Vol. 69, ed. E. Corbelli, D. Galli, and F. 

Palla (Florence: Soc. Ast. Italiana), p. 463 
Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, M., 

Pettini, M. 1999, ApJ, 519, 1 
Stpidel, C, C, Pettini, M., & Adelberger, K. L. 2000, ApJ, in press, 

fa,stro-ph/0008283| 
Stern, JJ., Bpmrad, H., Eisenhardt, P., Bunker, A. J., Dawson, S., 

Stanford, S. A., Elston R. 2000, ApJL, 533, 75 
Storrie-Lombardie, L. J., McMahon, R. G., Irwin, M. J., & Hazard, 

C. 1994, ApJ, 474, L13 
Weymann, R. J., Stern, D., Bunker, A., Spinrad, H., Chaffee, F. H., 

Thompson, R. I., & Storrie-Lombardi, L. J, 1998, A pi, 505,, 95 



Wood, K. & Loeb, A. 1999, ApJ, submitted, astro-ph/9911316 



